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This talk presents recent calculations in full QCD of the lowest three moments of generalized parton distributions 
and the insight they provide into the behavior of nucleon electromagnetic form factors, the origin of the nucleon 
spin, and the transverse structure of the nucleon. In addition, new exploratory calculations in the chiral regime 
of full QCD are discussed. 



1. INTRODUCTION 

Lattice field theory not only offers the prospect 
of precisely calculating the experimentally ob- 
servable properties of the nucleon from first prin- 
ciples, but also offers the deeper opportunity of 
obtaining insight into how QCD actually works 
in producing the rich and complex structure of 
hadrons. Beyond simply calculating numbers, we 
would like to answer basic questions of hadron 
structure. For example, how does the nucleon 
quark and gluon structure produce the observed 
scaling behavior of form factors? How does the 
total spin of the nucleon arise from the spin and 
orbital angular momentum of its quark and gluon 
constituents? What is the transverse, as well as 
longitudinal structure of the nucleon light-cone 
wave function? As the quark mass is continuously 
decreased from a world in which the pion mass is 
1 GeV to the physical world of light pions, how 
does the physics of the quark model and adiabatic 
flux tube potentials evolve into the physics of chi- 
ral symmetry breaking, where instantons, quark 
zero modes, and the associated pion cloud play 
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a dominant role? As discussed below, contempo- 
rary lattice calculation arc beginning to provide 
insight into these and other fundamental ques- 
tions in hadron structure. 

Because of asymptotic freedom, high energy 
lepton scattering provides precise measurements 
of matrix elements of the light-cone operator 

pA/2 

where n is a unit vector along the light-cone. 

The familiar quark distribution q(x) specifying 
the probability of finding a quark carrying a frac- 
tion x of the nucleon's momentum in the light 
cone frame is measured by the diagonal nucleon 
matrix element, (P\0(x)\P) = q(x). Expanding 
0{x) in local operators via the operator product 
expansion generates the tower of twist-two oper- 
ators, 

0W-^»}=^ 7 W B M2 ...i 5 M " } w (i) 

and the diagonal matrix element 
(P\0 { q ^^-^ } \P) specifies the (n - \) th mo- 
ment of the quark distribution J dx x n q(x). 
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The generalized parton distributions H(x, £, t) 
and E(x,£,t) |1I2I3I4| are measured by off- 
diagonal matrix elements of the light-cone opera- 
tor 

(P'\0(x)\P) = ((fi))H(x,£,t) 

+ l -^{{a av n a ))E{x^,t), (2) 

where A' 1 = P'f - P^, t = A 2 , £ = -n • A/2, 
and ((r)) = U{P')TU{P) for Dirac spinor U. Off- 
diagonal matrix elements of the tower of twist- 
two operators (P'\0^ 1 ' i2 '"' in ' i \P) yield moments 
of the generalized parton distributions, which in 
the special case of £ = 0, are 

J dxx^Hix^t) = A n<0 (t) 

J dxx n - 1 E{x,Q,t) = B n , {t), (3) 

where A n ^{t) and B n ^(t) are referred to as gen- 
eralized form factors (GFF's). Analogous ex- 
pressions in which the light-cone operator 0(x) 
and twist-two operators contain an additional 75 
measure the longitudinal spin density, Aq(x) and 
spin-dependent generalized parton distributions 
H(x, £, t) and E{x, £, t) with moments A n ^(t) and 
B n ,i(t). 

In this talk, I will discuss recent calcula- 
tions 5 6 7 8 of the generalized form factors 
-A(n=i,2,3),o(f) and A( n=li2 ,3),o(i) in full, un- 
quenched QCD in the currently computationally 
accessible domain that I refer to as the "heavy 
pion world" and discuss their physical signifi- 
cance. In addition, I will discuss initial efforts 
to explore the chiral regime in which the pion 
mass is sufficiently light that one can use chiral 
perturbation theory to extrapolate to the physi- 
cal pion mass. Although results in the heavy pion 
world cannot be directly compared with experi- 
ment, they nevertheless provide important insight 
into how QCD works, and provide the first step 
in the ultimate program of studying how hadronic 
physics evolves from the heavy pion world to our 
physical world. 



2. LATTICE CALCULATION 

The lowest three moments of spin-independent 
GPD's considered in this talk are 

(P'\O^\P) = ((r i ))A 10 (t) 
2m 

(p'\0^^\P)=P^(( 1 ^))A 2 Q(t) 

+ 7^P {AI1 (K 2}a »A Q S 20 (t) 
2m 

+ -A^A^CaoO), 
m 

+ ^-P^P^((a^ a ))A a B 3Q (t) 
2m 

+ A^A^<(7^>))A 3 2(t) 

+ ^-A^A^((a^ a ))A a B 32 (t), (4) 
2m 

Generalized form factors -A(n=i.2,3).oM ano - 
A( n =i,2,3),o(t) were calculated using the new 
method introduced in Ref. We consid- 

ered all the combinations of P and P' that 
could produce the same four-momentum trans- 
fer t — (P' — P) 2 , subject to the conditions 
that P = ^-(n x ,n v ,n z ) and P' = (0,0,0) or 
^-(—1,0,0). Using all these momentum combi- 
nations for a given t below 3.5 GeV, we calcu- 
lated all the H(4) cubic group lattice operators 
and index combinations producing the same con- 
tinuum GFF's, obtaining an overdetermined set 
of equations from which we extracted the most 
statistically accurate measurement of the desired 
GFF's the available lattice data can provide. As 
discussed in connection with Fig.0 the errors are 
substantially smaller than obtained by the com- 
mon practice of measuring a single operator with 
a single momentum combination. 

We calculated connected diagram contributions 
using approximately 200 SESAM 9 full QCD 
configurations with Wilson fcrmions at (3 = 5.6 
on 16 3 x 32 lattices. These calculations in the 
heavy pion world were performed at each of three 
quark masses, k— 0.1570, 0.1565, and 0.1560, cor- 
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Figure 1. 
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Electromagnetic form factor ratio 
The upper plot shows lattice 



results plotted with A=300 MeV. Experimental 
data taken from Ref. are shown below, and 
should be multiplied by the anomalous magnetic 
moment, K = 1.79, for comparison with the lattice 
results. 

responding to pion masses defined by ro of 744, 
831, and 897 MeV respectively. 

3. ELECTROMAGNETIC FORM 
FACTORS 

One of the early successes of perturbative QCD 
was the understanding of how the short range 
quark structure of a hadron governs the behav- 
ior of exclusive processes at large momentum 
transfer. However, whereas simple counting rules 
suggested that F 2 ~ Fi/Q 2 , experimental data 
from JLab shows that F 2 falls off much more 
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Figure 2. Generalized form factors A^ A {f) and 
i?2o d (£), with dipole fits denoted by dashed 
curves. 

slowly. Theoretically, it has recently been shown 
[TO] that the next to leading order light cone wave 
function yields F 2 ~ F 1 \og 2 (Q 2 /A 2 )/Q 2 , and the 
agreement between this prediction and the JLab 
data is shown in the lower portion of Fig. ^ 

Since the short range quark structure domi- 
nates this physics, it is reasonable to expect that 
omission of the pion cloud in the heavy pion world 
should not destroy the qualitative behavior. In- 
deed, our lattice results plotted in the top por- 
tion of Fig. \\\for the value A = 0.3 GeV yields 
excellent agreement with the Q 2 behavior of the 
experimental data. 

4. ORIGIN OF THE NUCLEON'S SPIN 

In the nonrelativistic quark model with three 
quarks in angular momentum zero single particle 
states, the total proton spin of 1/2 arises triv- 
ially from adding the spins of the three quarks. 
The so-called spin crises arose when deep inelas- 
tic scattering measurements of the lowest mo- 
ment of the spin-dependent structure function, 
AS = (1)a« + (1)a<j, indicated that only of the 
order of 30% of the nucleon spin arises from quark 
spins. Physically, in the heavy pion world where 
the quarks become less relativistic and are de- 
scribed by the quark model, it is reasonable to 
expect most of the nucleon spin to arise from the 
valence quark spin. As the quarks become lighter, 
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Figure 3. Example of the effectiveness of the 
overdetermined fit for C%f(t). The N=0 point 
uses three operators at a single external momen- 
tum combination to determine three form factors. 
The remaining points use six operators and N ex- 
ternal momentum combinations. 

one expects this fraction to decrease as the rela- 
tivistic quarks acquire more angular momentum 
and instanton effects, for example, remove helic- 
ity from the valence quarks and transfer it to the 
gluons and to quark- antiquark pairs. We would 
hope to observe this behavior clearly on the lat- 
tice. 

The total quark contribution to the nucleon 
spin is given by the extrapolation to t = of 
Al£ d (i) and B^ d (t) shown in Figure 2. Since 
d (t) is calculated directly at t = and B^it) 
is well fit by a constant that is measured to be 
nearly zero with small errors, the connected con- 
tribution to the angular momentum is measured 
to within a few percent. 

The calculation of the generalized form factors 
^20 (t), B2o{t), and C2o{t) provides an excellent 
example of the power of the new method intro- 
duced in Ref. 5 . Figure shows how adding ad- 
ditional operators and external momentum com- 
binations reduces the error in C2o{t) by a factor 
of 5 from the standard minimally determined case 
to our highly overdetermined case. 

Combined with the results of £ from reference 
|12|. we obtain the connected diagram contribu- 
tions to the decomposition of nucleon spin shown 



Figure 4. Fraction of the proton spin arising from 
the quark spin, AS, and from the quark orbital 
angular momentum, L q , as described in the text. 
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0.1570 


0.1565 


0.1560 


AS 


0.67± .04 


0.73 ± .03 


0.68 ± .02 


2L q 


0.06± .05 


-0.04 ± .04 


0.00 ± .03 


2J q 


0.73± .04 


0.69± .02 


0.68 ± .03 



Table 1 

Fraction of nucleon spin arising from quark spin, 
AS, quark orbital angular momentum, 2L q , and 
quark total angular momentum, 2J q . 



in Tab. and plotted in Fig. 21 Similar results 
have been obtained in Refs. |13I14| . To the extent 
that the disconnected diagrams do not change the 
qualitative behavior, we conclude that of the or- 
der of 70% of the spin of the nucleon arises from 
the quark spin and a negligible fraction arises 
from the quark orbital angular momentum in a 
heavy pion world where ~ 700 - 900 MeV. 
This behavior is just as expected from the ar- 
guments above. As lattice calculations approach 
the chiral limit, it will be interesting to fill in this 
graph and observe the quark spin contribution 
decrease to ~ 30% to agree with experiment. 



TRANSVERSE 
THE NUCLEON 



STRUCTURE OF 



In general, H (x, £, i) is complicated to inter- 
pret physically because it combines features of 
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Figure 5. Normalized generalized form factors 
A^ (t) and A££(t) forn=l (diamonds), n=2 (tri- 
angles) and n=3 (squares). 



Figure 6. Normalized generalized form factors 
A%* (t) and A^ (t) for n=l (diamonds), n=2 (tri- 
angles) and n=3 (squares). 



both parton distributions and form factors, and 
depends on three kinematical variables: the mo- 
mentum fraction x, the longitudinal component 
of the momentum transfer £, and the total mo- 
mentum transfer squared, t. In the particular 
case in which £ = 0, however, Burkardt JS] 
has shown that H(x,0,t), as well as its spin- 
dependent counterpart H(x,0,t), has a simple 
and revealing physical interpretation. 

It is useful to consider a mixed representation 
in which transverse coordinates are specified in 
coordinate space, the longitudinal coordinate is 
specified in momentum space and one uses light- 
cone coordinates for the longitudinal and time di- 
rections: x ± = (a;°±x 3 )/V2, p ± = {p°±p 3 )/V2. 
Using these variables, letting x denote the mo- 
mentum fraction and b± denote the transverse 
displacement (or impact parameter) of the light 



cone operator relative to the proton state, one 
may define an impact parameter dependent par- 
ton distribution 

q(x,b ± ) = {P + ,R_L= 0, A|G(a;,fe ± )|P + ,i? ± = 0, A) , 
where 

O q (x, b ± ) = J d ^-^\{- X —, b± ) 7 + q (?-, bl _). 

Burkardt shows that the generalized parton 
distribution H(x,0,t) is the Fourier transform of 
the impact parameter dependent parton distribu- 
tion, so that 

H{x,0,-A 2 ± ) = J d 2 b ± q{x,b ± )e lS± - K± 
A nfi (-Al) = J d 2 b ± Jdxx^ 1 

xq(x,b ± )e^ K \ (5) 
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tation, with the result ^B] 



q(x,b ± ) = (^) N - l J2J2 j 



3=1 



iS I 1 — y^ffj J 8 2 I a^rij J 5 (a; - x a ) 
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Figure 7. Transverse rms radius of the proton 
light cone wave function as a function of the aver- 
age quark momentum fraction, x av , for each mea- 
sured moment. 



where the second form follows from Eq. |J3}. Al- 
though one normally only expects a form factor to 
reduce to a Fourier transform of a density in the 
non-relativistic limit, Ref. ^H] shows that special 
features of the light cone frame also produce this 
simple result in relativistic field theory. Thus, if 
one could measure H(x,0,t) as a function of t, 
one could determine how the transverse distribu- 
tion of partons varies with momentum fraction 
x. 

Physically, we expect the transverse size of the 
nucleon to depend significantly on the momentum 
fraction x. Averaging q(x,b±) over all x, which 
produces Ai^(t) and thus corresponds to calcu- 
lating the form factor, the size is characterized by 
the transverse rms radius (r^_)5 = (x\ + x\)^ = 

|(r 2 )^. From the experimental electromag- 
netic form factor, the transverse rms charge ra- 
dius of the proton is 0.72 fm. As x — ► 1, the active 
parton carries all the momentum and the specta- 
tor partons give a negligible contribution. In this 
case the active parton represents the (transverse) 
center of momentum, and the distribution in im- 
pact parameter reduces to a delta function S(b±) 
with zero spatial extent. Indeed, explicit light 
cone wave functions |16I17| bear out this expec- 



where a denotes the index of the active parton, TV 
is the number of partons in the Fock state and the 
sum over c represents the sum over all additional 
quantum numbers characterizing the Fock state. 
Here, one explicitly observes lim x ^i q(x, b±) oc 
5(b±). Since H(x, 0, t) is the Fourier transform of 
the transverse distribution, the slope in —t = 
at the origin measures the rms transverse radius. 
As a result, we expect the substantial change in 
transverse size with x to be reflected in an equally 
significant change in slope with x. In particu- 
lar, as x — > 1, the slope should approach zero. 
Hence, when we calculate moments of H(x,0,t), 
the higher the power of x, the more strongly large 
x is weighted and the smaller the slope should be- 
come. Therefore, this argument makes the quali- 
tative prediction that the slope of the generalized 
form factors A ni o(i) and A n fi{t) should decrease 
with increasing n, and we expect that this effect 
should be strong enough to be clearly visible in 
lattice calculations of these form factors. 

Figures and show the generalized form fac- 
tors A U} o(t) and A n ^(t) for the lowest three mo- 
ments, n = 1, 2, and 3. The form factors have 
been normalized to unity at t = to make the de- 
pendence of the shape on n more obvious. Note 
that ^4i,o, ^3,0j and ^2,0 depend on the difference 
between the quark and antiquark distributions 
whereas Ai t o,Asfi, and ^2,0 depend on the sum. 
Hence only comparisons between moments differ- 
ing by n — 2 compare the same physical quantity 
with different weighting in x. To facilitate de- 
termination of the slope of the form factors and 
to guide the eye, the data have been fit using a 
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dipole form factor 

-dipole _ a / fi s 

A nfl ~ (I- 4 )2 ' ^ > 

The solid line denotes the least-squares fit and the 
shaded error band shows the error in the slope 
Arrid given by the fit. Although the dipole fit is 
purely phenomenological, we note that it is sta- 
tistically consistent with the lattice data. 

The top panel in Fig. [S] shows the flavor non- 
singlet case A u — A d , for which the connected 
diagrams we have calculated yield the complete 
answer. It is calculated at the heaviest quark 
mass we have considered, corresponding to m, 
= 0.87 GeV. Note that the form factors are sta- 
tistically very well separated, and differ dramati- 
cally for the three moments. Indeed, as discussed 
more quantitatively below, the slope at the ori- 
gin decreases by more than a factor of 2 between 
n = 1 and n = 3, indicating that the transverse 
size decreases by more than a factor of 2. The 
top panel of Fig. shows analogous results for 
lighter quarks, — 0.74 GeV, where we ob- 
serve the same qualitative behavior but slightly 
weaker dependence on the moment. The second 
panel of Fig. shows the flavor singlet combina- 
tion A u + A d , for which we have had to omit the 
disconnected diagram because of its significantly 
greater computational cost. Comparing this fig- 
ure with the top panel calculated at the same 
quark mass, we note that while the connected 
contributions to A u ±A d are qualitatively similar, 
there is significant quark flavor dependence that 
can be used to explore the nucleon wave func- 
tion. The bottom panel of Fig shows the spin- 
dependent flavor non-singlet form factors A u — A d 
at the heaviest quark mass. Comparison with the 
top of Fig. [S] displays the difference between the 
spin averaged and spin dependent densities. We 
observe a striking difference, in that the change 
between the n = 1 and n = 3 form factors for 
q{x, — q(x, b±)i is roughly 6 times smaller 
than for \{q{x, b±)^ + q(x, b±)i). 

Finally, it is useful to use the slope of the form 
factors at t = to determine the transverse rms 
radius, 

2 {n) J^b^lfdxx^qjx^^) 
K ±> Jd 2 b ± Jdxx n - 1 q(x,b ± ) ' 1 ' 



Transverse rms radii calculated in this way for 
the first three moments are plotted in Fig. for 
m x = 0.870 GeV. To set the scale, the transverse 
charge radius at this mass is (r^) c h ar g e = 0.48 fm, 
which is two-thirds the experimental transverse 
size 0.72 fm, reflecting the effect of the absence 
of a significant pion cloud. The nonsinglet trans- 
verse size (r^) u _d = 0.38 fm is slightly smaller 
than the rms charge radius, but drops 62% to 
0.14 fm for n=3. The singlet size (r^) u+( j is 0.46 
fm, and drops 43% to 0.27 for n=3. This is a 
truly dramatic change in rms radius arising from 
changing the weighting by x 2 . 

To display the change in transverse size with 
x in Fig. the rms radii for each moment are 
plotted at the average value of x corresponding to 
that moment. If one neglects the fact that the an- 
tiquark distribution contributes to even and odd 
moments with opposite sign, the mean value of x 
in the distributions q{x), xq{x), and x 2 q{x) are 
determined directly from the moments of struc- 
ture functions measured on the lattice. Applying 
a small correction for the effect of the alternating 
antiquark contributions yields the mean values of 
x for each moment plotted in the Figure. The x 
dependence shown in this Figure is quite striking, 
with the nonsinglet transverse size dropping 62% 
as the mean value of x increases from 0.2 to 0.4, 
and going to zero when x reaches 1.0. 

6. EXPLORATORY CALCULATION IN 
THE CHIRAL REGIME 

Full QCD calculations with light quark masses 
are notoriously expensive, so significant compro- 
mises are required to begin to explore the chi- 
ral regime. Our initial exploration of this regime 
is a hybrid calculation using MILC configura- 
tions []3| with staggered sea quarks and domain 
wall valence quarks. The MILC configurations 
on a 20 3 x 64 lattice use strange quark masses 
am s = 0.05 and light quark masses am u +d = 0.01 
and 0.05 with the Asqtad action corresponding 
to lattice spacing 0.13 fm. HYP-smearing [2"U| 
with ol\ — 0.75, a2 = 0.6, and 03 = 0.3 was used 
to reduce the effect of dislocations. Chiral va- 
lence quarks were calculated using domain wall 
fcrmions with L$= 16 and M — 1.7. The lat- 
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Figure 8. Hybrid chiral calculation of (x) l 
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tice size of 2.6 fm can sustain pions as light as 
300 MeV, and our initial calculations were carried 
out for pion masses of approximately 343 and 635 
MeV. 

In the long term, our plan is to attain high 
statistics on several lattice volumes and two cou- 
pling constants for a range of quark masses. Hy- 
brid partially quenched chiral perturbation the- 
ory will be used to correct for the inconsistency 
between chiral valence quarks and staggered sea 
quarks, and perturbative renormalization is ex- 
pected to be adequate because of the improved 
convergence arising from HYP-smearing. 

However, the initial explorations shown in 
Figs. |H1 and [5] are based on roughly 100 configu- 
rations, only have tree-level renormalization, and 
are not corrected using chiral perturbation the- 
ory. Figure |H1 shows the well-known case of the 
momentum fraction, (x), for which the chiral ex- 
trapolation formula |23 is roughly constant in the 
heavy pion regime of the three SESAM points, 
and then decreases sharply in the region of m\ ~ 
0.2 GeV 2 by 50% to agree with the experimental 
point. The two MILC points, which unlike the 
renormalized SESAM results must still be renor- 
malized, do not yet give any indication of decreas- 
ing in the chiral regime, and it is an open ques- 
tion as to whether the lowest point is subject to 
large finite volume corrections. The axial charge 
shown in the upper panel of Figgis also presently 
unsatisfactory at the light quark point. In this 
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Fi gure 9. Hybrid chiral calculation of (1)a«— Ad 
and {x)au- Ad- 



case, it is well known that finite volume effects 
produce large discrepancies, so calculations on a 
larger lattice are clearly needed in this case. Only 
for the case of the first moment of the spin dis- 
tribution, (x) Au-Ad, is the qualitative behavior 
roughly consistent with experiment. 

Whereas these initial results are still far 
from being quantitatively controlled, they clearly 
demonstrate feasibility of hybrid calculations in 
the chiral regime. Improvement of statistics, care- 
ful study of finite volume effects, and calcula- 
tion of hybrid chiral perturbation corrections and 
renormalization factors offer the potential for a 
first glimpse at hadron structure in the chiral 
regime. 
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7. SUMMARY AND OUTLOOK 

In the heavy pion world presently accessible to 
unquenched lattice QCD, we have calculated the 
lowest 3 generalized form factors A n fl and A n $ 
to — t — 3 GeV 2 and shown that they provide 
insight into several important aspects of hadron 
structure. We note that in cases that are com- 
parable, our results are consistent with the cal- 
culations of the lowest two moments in Ref. |13| . 
The overdetermined method for measuring gener- 
alized form factors produces good statistics up to 
3 GeV 2 , and enables meaningful study of electro- 
magnetic form factors, calculation of the origin 
of the nucleon spin, and study of the transverse 
structure of light cone wave functions. A particu- 
larly striking result is the dramatic 62% decrease 
in the transverse size (^5_) u _d between the first 
and third moment. We also observed clear de- 
pendence of the transverse distribution on flavor 
and spin and have shown that the commonly used 
factorization Ansatz H(x, 0, t) — Q(x)F(t) is fun- 
damentally wrong. 

The most immediate challenges are to extend 
these calculations to the chiral regime of realistic 
quark masses, and to extend techniques for evalu- 
ating disconnected diagrams [22 to these observ- 
ables. When precise, controlled extrapolations to 
the physical pion mass are finally achieved, they 
will play a special role in our understanding of 
hadron structure. 
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